Dynamics and phase evolution of Bose-Einstein condensates in one-dimensional 

optical lattices 



O. Morsch*, M. Cristiani, J.H. Miiller, D. Ciampini, and E. Arimondo 

INFM, Dipartimento di Fisica E. Fermi, Universita di Pisa, Via Buonarroti 2, 1-56127 Pisa, Italy 

* e-mail: morsch@df.unipi.it 
(February 1, 2008) 



O , 
O . 
(N 

^■ 
CD 

(N 



> 

(N 
O 
0^ 
O 
(N 
O 



X3 
O 

o 



X 



We report experimental results on the dynamics and phase 
evolution of Bose-Einstein condensates in ID optical lattices. 
The dynamical behaviour is studied by adiabatically loading 
the condensate into the lattice and subsequently switching 
off the magnetic trap. In this case, the condensate is free 
to expand inside the periodic structure of the optical lattice. 
The phase evolution of the condensate, on the other hand, 
can be studied by non-adiabatically switching on the periodic 
potential. We observe decays and revivals of the interference 
pattern after a time-of-flight. 

PACS number(s): 03.75.Fi,32.80.Pj 



I. INTRODUCTION 

The study of Bose-Einstein condensates (BECs) in pe- 
riodic potentials has recently seen some major advances 
both on the theoretical [1-4] and, most notably, the ex- 
perimental side [5,6]. Given the rapid progress made 
in optical lattices research with ultra-cold atoms in the 
1990s, it was only a matter of time before the BEG com- 
munity would extend those studies to the domain of Bose- 
condensed atoms. 

So far, experiments with Bose-Einstein condensates 
have focused mainly on two areas. On the one hand, 
there have been 'conventional' investigations of the dy- 
namics of condensates in stationary, accelerated and 
pulsed optical lattices [7,8]. On the other hand, the very 
efficient loading of ultra-cold atoms into a periodic struc- 
ture afforded by BEGs has also made possible the obser- 
vation of quantum effects such as number squeezing [5] 
and the Mott-insulator transition [6] . 

In previous experiments, we have already addressed 
the issue of a BEG in an accelerated optical lattice and 
associated phenomena such as Bloch oscillations and 
Landau-Zener tunneling [7,9]. In the present work, how- 
ever, we shall concentrate on the properties of BEGs in 
stationary periodic potentials. We have investigated both 
the dynamics and the phase evolution of BEGs in such 
potentials. 

This paper is arranged as follows. After briefly pre- 
senting our experimental setup in section II, we discuss 
the results of our experiments on the free expansion of 
the BEG in the lattice (section III) . The phase evolution 
of the condensate is the subject of section IV, which is 
followed by conclusions and an outlook on further exper- 



iments (section V). 



II. EXPERIMENTAL SETUP 

Our experimental setup is described in detail else- 
where [7,9]. Briefly, after creating BEGs of iVo = 1 — 
2 X 10'' *^Rb atoms in a triaxial time-orbiting potential 
trap, we adiabatically lower the trap frequency V^^ap to 
the desired value and then superimpose onto the mag- 
netic trap an optical lattice (either along the vertical or 
the horizontal trap axis) created by two linearly polar- 
ized Gaussian laser beams intersecting at a half-angle 9 
and detuned by « 30 GHz above the rubidium resonance 
line. The periodic potential U{z) — Uq srn^ {tt z / d) thus 
created [10] has a lattice spacing d = 7r/[fcsin(6'/2)], and 
the depth Uq of the potential (measured in lattice re- 
coil energies Ej-ec = h^TT^ /2md'^) can be varied between 
and ~ 20 Erec by adjusting the laser intensity using 
acousto-optic modulators. In our setup, we could real- 
ize a horizontal lattice with d = 1.56 /im and a vertical 
lattice with d = 1.2 /im. The number of lattice sites occu- 
pied by the condensate lay between 10 and 15, depending 
on the trap frequency and the lattice spacing. When us- 
ing the vertical lattice with the magnetic trap switched 
off (in the expansion experiments, see below), we accel- 
erated the lattice downwards by chirping the frequency 
difference between the lattice beams in order to offset 
gravity in the rest-frame of the lattice. 



III. DYNAMICS: FREE EXPANSION INSIDE A 
ID LATTICE 



We studied the dynamics of the condensate in a sta- 
tionary lattice by observing the free expansion of a con- 
densate adiabatically loaded into the lattice. In order 
to load the condensate into the optical lattice, after adi- 
abatically reducing the mean magnetic trap frequency 
i^trap to the desired value we linearly increased the lat- 
tice depth from to Uq in a time tramp = 150 ms. Since 
typically the chemical potentials /xq of our condensates 
in traps with frequencies Vtrap/^iT « 20 — 70 Hz lie be- 
tween 50 Hz and 200 Hz, the adiabaticity condition [11] 
tramp > h / {Iq was Satisfied even for small chemical po- 
tentials. After that, only the magnetic trap was switched 
off and the condensate was imaged after a variable time 
tf^xp of free expansion inside the optical lattice. The effect 
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of a deep lattice on the free expansion is clearly evident 
in the condensate images of Fig. 1 (a) and in the mea- 
surements of the condensate aspect ratio p\\l ps_ shown in 
Fig. 1 (b) (p||,_L denoting the half-widths of Gaussian 
fits to the density profile in the lattice direction and per- 
pendicular to it, respectively). As in the presence of the 
lattice the condensate does not expand at all in the lattice 
direction, the aspect ratio drops off sharply with increas- 
ing time-of-flight as the BEG expands in the (unbound) 
perpendicular direction. The lack of expansion in the 
lattice direction reflects the fact that the condensate has 
effectively been split up into several smaller condensates 
confined in the individual lattice wells. In the perpendic- 
ular direction we have observed an enhanced expansion 
when the lattice is present, which can be explained by 
the increase in the chemical potential when the lattice is 
ramped up (as calculated by Pedri et al. [3]). 

We modeled the expansion of the condensate in the 
presence of the ID lattice using the equations derived 
in [12,13] with a slight modification along the linos of 
Ref. [4]: In the ansatz leading to the differential equa- 
tion of [12] for the scaling factor A||(i = t^xp) ~ P\\{t = 
texp) /P|| (i = 0) along the lattice direction we replaced the 
atomic mass m by the effective mass m*{Uo), introduced 
in [4] in analogy with a solid state physics approach, as 
derived from a band structure calculation for a periodic 
potential of depth Uq. As can be seen in Fig. 2, taking 
into account the variation of the chemical potential with 
Uo this model reproduces well our experimental data for 
the perpendicular expansion of the condensate (the the- 
oretical plots are corrected for the 5 /xm resolution of our 
imaging system). 

For the lattice direction, this approach gives the cor- 
rect result for deep lattices (Uo ^ Erec), for which intu- 
itively one expects the condensate to be broken up and 
the individual parts to be essentially confined to the lat- 
tice wells, suppressing the expansion along the lattice 
direction. The qualitative behaviour in the intermediate 
regime (0 < f/o < 5 Erec) is also reasonably well repro- 
duced. Wo find, however, that the experimental expan- 
sion along the lattice direction is considerably less than 
the theoretical prediction in the intermediate regime. We 
have checked that taking into account the finite momen- 
tum spread of the condensate when calculating the ef- 
fective mass only leads to a correction on the percent 
level and thus cannot explain the deviation of our exper- 
imental data from the numerical calculations neglecting 
mean field corrections. This might indicate that there 
are effects such as self-trapping [2] due to the mean-field 
interaction that further reduce the expansion in the lat- 
tice direction. Another hint in this direction is given by 
the fact that when we vary the number of atoms N in the 
condensate, for large N the width in the lattice direction 
of the expanded condensate actually starts to decrease 
rather than increase. Again, this might be explained by 
mean-field effects which become important in the dynam- 
ics of the condensate expansion when N gets large, with 



non- linear effects like self-trapping reducing the observed 

width p|| . 

As described in detail in [14], we could also directly 
deduce from our data the dependence of the chemical 
potential on Uq from the initial perpendicular size p±_ {t = 
0) of the condensate in the presence of the lattice inferred 
from the size p±{t = texp) measured after an expansion 
time texp and the ratio X±{t = texp)/X±{t = 0) of the 
scale factor A (t) . The chemical potential thus measured 
agrees well with the prediction of [3] . 

Finally, we note here that in the limit of large lattice 
depths, the adiabatic loading of a BEG into and optical 
lattices effectively realizes an adiabatic transformation 
between a 3D condensate and an array of 2D conden- 
sates. The condition fiso < TiLOiat of Ref. [15] (where 
^lat = '^{Erec/K}\/Uo / Erec is the harmonic approxima- 
tion for the oscillation frequency in a lattice well) for 
the condensates in each well to be in the 2D limit is al- 
ways satisfied for the small number of atoms in a single 
well (« 10'^) present in our experiment. For an array 
of 2D condensates obtained by creating the condensate 
in the combined potential of the harmonic trap and the 
lattice. Burger et al. [16] have shown that in the case of 
their cigar-shaped condensate (with the long axis along 
the lattice direction), the transition temperature T^^ in 
the presence of the lattice is significantly lower than T^^ 
in the 3D case (i.e. in the magnetic trap without the 
lattice). Galculating the critical temperature T^-° along 
the same lines for our system, we find that T^^ « T^^ 
due to the larger number of atoms per lattice site in our 
geometry, and hence we expect no significant change in 
the condensate fraction in the presence of the lattice. In 
fact, experimentally we even find a consistently larger 
condensate fraction after ramping up the lattice, as seen 
in Fig. 3. This result indicates that, with an appropriate 
choice of parameters, a ID optical lattice could be used 
to investigate adiabatic transformations between 3D and 
2D condensates which could, e.g., be exploited to create 
condensates from thermal clouds by changing the dimen- 
sionality of the system, as recently proposed by Olshani 
and Weiss [19], similarly to the change in the shape of 
the potentials in Refs [17,18] for other geometries. 

IV. PHASE EVOLUTION: BEHAVIOUR OF THE 
INTERFERENCE PATTERN 

In order to study the phase evolution of a condensate in 
an optical lattice, we performed a non- adiabatic loading 
procedure, meaning that the condensate initially was not 
in the ground state of the combined potential of the mag- 
netic trap and the optical lattice. In a typical experiment, 
the optical lattice was ramped up in Tramp « 1 — 5 ms, af- 
ter which the potential was kept at its maximum value Uq 
for a holding time thoid- At the end of the holding time, 
the lattice was accelerated in 1 ms to one (lattice) recoil 
velocity by chirping the frequency difference between the 
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lattice beams. Immediately after that, both the lattice 
and the magnetic trap were switched off. After a time-of- 
flight of 20 — 22 ms, the expanded condensate was imaged 
using a resonant probe flash. Figure 4 shows typical ab- 
sorption images and integrated lines profiles obtained in 
this manner for different holding times. For short times, 
a clean double-peak structure is visible, as expected from 
the interference between the condensates expanding from 
the individual lattice wells (with a tt phase difference be- 
tween them due to the final acceleration). During the 
first few milliseconds, this pattern evolves into a more 
complicated structure featuring several additional peaks, 
and finally washes out completely, resulting in a single 
Gaussian-shaped lump. For long waiting times, the two- 
peaked structure reappears. 

Intuitively, these results can interpreted as follows. Ini- 
tially, the condensate is broken up into several pieces 
when the optical lattice is abruptly switched on, and 
the condensate fragments are locally compressed. This 
leads to a different (local) chemical potential and hence 
a different phase evolution at the individual lattice sites. 
As a result, the interference pattern after a time-of- flight 
evolves from a clean two-peaked structure (all condensate 
fragments in phase) into a broad Gaussian distribution, 
reflecting the fact that the individual fragments have ac- 
cumulated different phases. For long holding times, a 
common phase is re-established through tunneling be- 
tween adjacent lattice sites and dissipation of energy into 
higher-lying modes (more detailed studies of these phe- 
nomena are under way). 

In order to quantify the degree of coherence between 
the lattice sites, we characterized the interference pat- 
tern (integrated perpendicular to the lattice direction) 
through a visibility ^ calculated as 

^ hpeak ^middle \ 

S — T TT 5 \^) 

•^peak ' '''middle 

where hpeak is the mean value of the absorption image 
at the position of the two peaks, and hmiddie is the value 
of the absorption image mid- way between the two peaks 
(averaged over a range of 1/5 of the peak separation) 
as illustrated in Fig. 4 (a). In a typical experiment, ^ 
initially rapidly decreases from a value of « 0.6 — 0.9 
to roughly within f« 5 — 20 ms, depending on the trap 
frequency Virap (sec Fig. 5). If the lattice depth is kept 
constant throughout, for trap frequencies Vtrap > 30 Hz 
the visibility then typically rises again and begins to fluc- 
tuate in an apparently random manner for a few tens of 
milliseconds. Finally, these fluctuations die out and ^ 
stabilizes at a value close to the initial visibility. 

For weak magnetic traps with Vtrap < 30 Hz, after 
the initial decay in visibility ^ stays close to for up to 
200 ms. Under these conditions, we could induce a faster 
revival of the interference pattern by lowering the opti- 
cal lattice depth and hence increasing the tunneling rate 
between adjacent lattice wells. Fig. 5 shows the initial 
decay for a lattice depth of 12 Erec and the subsequent 
revival when the lattice is lowered to w 4£^rec- 



The interpretation of the phase evolution is compli- 
cated by the fact that when the lattice is switched on 
non-adiabatically. radial oscillations of the condensate 
are also excited. These oscillations damp out after a few 
cycles. Also, the temperature of the condensate initially 
rises and levels out at a higher value. In our experi- 
ments, T/Tc (where Tc is the critical temperature for 
condensation) typically rises from 0.7 to « 0.85. This 
re-thermalization, which takes places on a timescale that 
is comparable both to the damping time of the radial 
oscillations and the revival time of the visibility at fixed 
lattice depth, may take place through a re-distribution of 
energy into higher-lying modes of the condensate. Theo- 
retical studies of these processes are currently being car- 
ried out. 



V. CONCLUSIONS AND OUTLOOK 

In a stationary optical lattice, both dynamical and 
phase properties of a BEG can be studied. By adiabat- 
ically loading the BEC into the lattice, we have studied 
its free expansion after switching off the magnetic trap. 
From the expansion perpendicular to the lattice direc- 
tion we were able to deduce the chemical potential of 
the BEC in the lattice as a function of the lattice depth. 
The expansion in the lattice direction for small lattice 
depths still needs further investigation as it is not well 
reproduced by a simple treatment involving the effective 
mass. Non-linear effects might come into play, possibly 
reducing the expansion through self-trapping. 

The phase evolution of a condensate non-adiabatically 
loaded into a lattice is another intriguing area of study 
that deserves closer investigation. We have observed an 
initial decay of the interference pattern followed, for long 
holding times, by an almost complete revival. The role 
of tunneling in this revival is evident from the fact that 
when the lattice depth is lowered after the initial dephas- 
ing period, a revival in the visibility is induced. The in- 
terplay between tunneling, re-thermalization and damp- 
ing of the radial oscillations excited by the non-adiabatic 
switch-on of the lattice needs to be looked at more closely 
in future investigations. In particular, the role of the non- 
condensed fraction needs to be elucidated. 
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FIG. 1. Effect of the one-dimensional optical lattice on the 
condensate expansion. In (a), the difference in condensate 

width along the lattice direction after ~ 23 ms of free expan- 
sion with (solid line) and without the lattice (dashed line) is 
clearly visible. The evolution of the aspect ratio of the con- 
densate in a lattice with Uo = 20i?rec is shown in (b). When 
the lattice is present (filled symbols), the condensate does not 
expand in the lattice direction and hence the aspect ratio de- 
creases. Without the lattice (open symbols), for the same 
magnetic trap the aspect ratio increases. 
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FIG. 2. Dependence on the lattice depth of the condensate 
widths after 22 ms time-of-flight. The initial size (including 
the 5 /um resolution of our imaging system) of the condensate 
is ~ 10 /ini. In the lattice direction, for lattice depths > 2 Erec 
the condensate docs not expand anymore. This early cut-off is 
not predicted by our simple scaling approach using the effec- 
tive mass (dashed line). The expansion in the perpendicular 
direction, however, is reproduced well by this model (solid 
line) when the increase in chemical potential with Uo is taken 
into account (see text). In this experiment the mean trap 
frequency was 25 Hz, with i'±/i'\\ = V^/l. 
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FIG. 3. The condensate fraction Nq/N with (filled sym- 
bols) and without (open symbols) optical lattice as a func- 
tion of the final rf-cut in the evaporation. When the lattice 
is present, the condensate fraction is comparable to or even 
slightly larger than without the lattice. In this experiment, 
Uo m 15 Erec and the mean trap firequency was 26 Hz. 
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FIG. 4. Evolution of the interference pattern of a con- 
densate released from an optical lattice after non-adiabatic 

loading (absorption images and integrated profiles). The dis- 
tinct two-peaked structure visible immediately after loading 
(a) washes out after a few milliseconds (b). For long hold- 
ing times in the lattice, the initial structure reappears (c). In 
this experiment, the mean trap frequency was 26.7 Hz and the 
lattice depth Uo ~ 15 Erec- The holding times in the lattice 
from (a) to (c) were 1,22 and 300 ms, respectively. In (a), the 
definition of the visibility is illustrated. 
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FIG. 5. Decay and revival of the interference pattern. Af- 
ter the initial dcpliasing in a lattice with Uo ~ 12 Erec (filled 
symbols), the lattice depth was lowered to ~ 4 Erec (open 
symbols), resulting in a faster revival of the visibility. The 
solid and dashed lines are exponential fits with time con- 
stants Tdephase ~ 18 ms and Trephase ~ 50 ms, respectively. 
The dash-dotted line represents the visibility observed for a 
condensate adiabatically loaded into the lattice with a ramp 
time Tramp = 150 ms. 
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